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ABSTRACT 


We discuss models for hot gas in the broad line region of active galactic 
nuclei. The results of the two phase equilibrium models for confinement of 
broad line clouds by Compton heated gas as described by Krolik, McKee and 
Tarter (1982) are used to show that high luminosity Quasars are expected to 
show Fe XXVI La line absorption which will be observable with spectrometers 
such as those planned for the future X-ray spectroscopy experiments. Two 
phase equilibrium models also predict that the gas in the broad line clouds 
and the confining medium may be Compton thick. We show that the combined 
effects of Comptonization and photoabsorption can suppress both the broad 
emission lines and X-rays in the Einstein and HEAO-1 energy bands. The 
observed properties of such Compton thick active galaxies are expected to be 
similar to those of Seyfert 2 nuclei. The implications for Dolarization and 
variability are also discussed. 


Subject Headings: galaxies: nuclei -- quasars --radiation mechanisms 



I. INTRODUCTION 


Several recent theoretical and observational developments suggest that 
hot (T > 10 6 K) gas may be important in the broad emission line reqion of 
quasars and active galactic nuclei (AGN). Theoretical arguments concerning 
the thermal instability of photoionized gases (McCray 1979; Krolik, McKee and 
Tarter, 1982), together with the need for pressure confinement of broad line 
clouds (e.g. Mathews 1974) predict that the clouds will be bathed in a hot 
confining medium. If this is so, then the hot intercloud qas (ICM) may 
produce observable X-ray absorption and emission lines from ions such as 
hydrogen-like iron, in addition to continuum absorption edges. 

Hot gas may also have broad implications in determining the 
distinguishing features of the various different classes of AGN. If the qas 
is present in slightly larger quantities than predicted by thermal instability 
(two phase) models then Comptonization can suppress the highest energy 
continuum from the central source. Such Comptonized spectra will not have an 
unstable effect on gas in the broad line region (the broad emission line 
clouds may not form), and the X-ray continuum may also be photoabsorbed in the 
Comptonizing region. Such Compton thick objects may resemble Seyfert 2's in 
their continuum properties. 

In this paper we consider the spectral consequences of hot gas in the 
broad emission line region of AGN. We do this both in the context of the two 
phase models proposed for objects which show broad lines, and in the more 
general case when the continuum source is embedded in Compton thick 
material. In doing so, we employ standard "nebular" approximations for the 
transfer of continuum photons when the gas is Compton thin (see, e.g. Kallman 
and McCray 1982; Kallman 1984). Gas flows and dynamical models for broad line 
regions are considered only in so far as they affect the gas distribution in 
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the broad line region. Narrow emission lines and radio frequency properties 
of AGN are not considered. Compton thin models are presented in $11. The 
Compton thick case is considered in §111. §IV provides a discussion and 
summary of our results. 

II. SPECTRAL FORMATION IN THE QUASAR ICM 


a) The Two Phase Model 

In this section we explore the potentially observable proDerties of the 
ICM in objects which are observed to have broad optical and UV emission lines, 
i.e. quasars and Seyfert 1 galaxies. In such objects the correspondence 
between the ICM and the properties of the broad line emittinq clouds has been 
discussed in detail by Krolik, McKee, and Tarter (1982; hereafter KMT). We 
present here a summary and confirmation of their results as a background for 
our later discussion. 

The connection between the broad emission line clouds and the ICM is due 
to the fact that both are heated and photoionized primarily by continuum 
photons from a compact source at the center of the broad line region. The qas 
is assumed to exist in thermal equlibrium at constant pressure, P. The gas 
temperature and degree of ionization depend only on the shape of ionizinq 
continuum spectrum and on a parameter proportional to the ratio of continuum 
radiation pressure to gas pressure, the ionization parameter H. In terms of 
the incident flux, F (erg sec -1 cm -2 ), the ionization parameter is defined as 
2 = F/cP. At large values of this parameter the gas is fully ionized and the 
temperature is determined by Compton and inverse-Compton scattering. The 
temperature in this case, denoted Tjq, depends on the shpae of the continuum 
spectrum and on the high and low energy boundaries of the spectrum. For 
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example, if the spectral shape is a power law with (energy) index less than 
one T IC is sensitive to the high energy spectral cutoff. If the spectral 
index is greater than one, Tic de P ends on the cutoff at low energies. At 
small H the gas is recombined, and the temperature is approximately 10 4 K, 
owing to the strong cooling by atomic transitions of hydroqen. Hell, and C IV 
ions. At intermediate S the thermal equilibrium condition may allow both 
Comptonized and 10 4 K temperatures to occur at the same value of S. In this 
two phase region of parameter space the net cooling function of the gas 
satisfies Field's (1965) criterion for thermal instability. The thermally 
unstable zone is bounded by the critical ionization parameters, denoted 

and s* (3£ > 2*). The value of E* is fixed by the atomic coolinq properties 
of the gas, =■* a 10. The value of 2* depends on T^: S£ a T^ 372 . The two 
phase zone is evanescent (i.e. 2£ < 2*) for T JC < 10 7 K. 

The two phase models are successful in accounting for the fact that 
quasar broad emission line clouds are observed to exist over a narrow range of 
ionization parameter. Comparison of observed broad line ratios with cloud 
photoionization models reveals that clouds exist in the range of ionization 
parameters 0.3 < S < 2 (Kwan and Krolik 1981) for an ensemble of AGN with 
luminosities in the range 10 44 erg s -1 - 10 47 erg s"*. This range is rouqhly 
congruent with the range of ionization parameters for which clouds and ICM can 
coexist via the two-phase models. It is well known from measurements of the 
line fluxes and sizes of the BLR that the filling factor of the dense clouds 
is «1. Therefore there is not enough gas to fill the BLR at 10 4 K. . The 
rates of cloud evaporation and free expansion are sufficient to produce the 
ICM even if the clouds were somehow initially produced in the absence of such 
hot gas. 
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The simple two phase model is subject to modification by a number of 
other processes which may be occuring in AGN. Sources of heating other than 
photoionization and Compton scattering, such as suprathermal particles or 
radio frequency heating, will increase the ICM temperature above Tj C and 
extend the range of thermally unstable ionization parameters. Advection of 
material into the broad line region by gas flowing with a timescale comparable 
to the cooling and heating timescale can act either to enhance the instability 
(if the advected material is hot) or to partially suppress it (if the advected 
material is cool). Evaporation of the clouds will tend to raise the cloud 
pressure over that of the ICM. Radiation pressure due to trapped line photons 
can disrupt clouds preferentially at large ionization parameters, thus 
effectively reducing (see also Weymann 1976). Although all these processes 
are potentially important in determining the ICM properties, a detailed 
accounting of them depends sensitively on many poorly known properties of the 
broad line clouds (gas flow, spectrum of cloud sizes, relativistic particle 
flux, etc.). In what follows we adopt the point of view that the simple two 
phase model provides an approximate average description of the ICM, and we use 
this for purposes of quantitative estimates of the ICM properties. We also 
mention the potential qualitative modification to our results that will arise 
from the processes listed in this paragraph. 

b) Two Phase Models 

The existence of observable spectral features formed in the ICM depends 

on the degree of ionization of the ICM, and on the gas column density. The 

7 ft 

ion abundances depend on the gas temperature, T « T IC * 10 -10 K, and on the 
shape of the ionizing continuum spectrum. The column density depends on the 
values of the critical ionization parameters, and E*, which also depend on 
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the continuum shape, and on the distribution of gas pressure as a function of 
position. For a variety of plausible choices for continuum shape and pressure 
distributions the two phase model predicts the existence of features which 
will be observable with future X-ray spectroscopy experiments. 

The value of the ICM temperature, Tj£, depends on the detailed shape of 
the continuum from the central source. Observed AGN spectra typically have a 
power law spectral index a « 0.7 (in enerqy) between 1 keV and several hundred 
keV (Mushotzky et al . 1980; Worrall et al. 1981; Primini et al. 1980; 
Rothschild et al. 1982) and a « 1.2 between wavelength of 10y and 1000 A 
(Malkan and Sargent 1980; Glass 1982, Glass et al. 1982). The spectra between 
0.1 eV and 500 keV can be adequately represented by 


f e = 


A e -1 * 2 e _£/e o + B e -0 * 7 


( 1 ) 


where e is the photon energy in eV and e Q = 200 eV. The spectral softness 
parameter A/B has a value in the range from 0.1 (eV) 1 / 2 to 1000 (eV) 1 / 2 
(Zamorani et al . 1981; Reichert et al. 1982). This ratio is related to the 
conventionally used optical to X-ray spectral index, o , according 

U A 

to a «v = 1°9 A/B + 0.698. For A/B = O.l(eV) 1 / 2 , a « 0.3; for A/B = 

OX Z.bUo ox 

1000 (eV) 1 / 2 , a - 1.8. 

0 X 

Figure 1 shows the range of E£, s*, and T^ as a function of the spectral 
softness parameter A/B. These results were calculated using the 
photoionization program described by Kallman and McCray (1982) and Kallman 
(1984) using the element abundances of Kwan and Krolik (1982). The incident 
spectrum is assumed to extend from 0.1 eV to 511 KeV. Compton heating and 
cooling are treated using the fully relativistic cross section and electron 
velocity distribution function. The range of broad line cloud ionization 
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parameters inferred from the ensemble of observed quasar line ratios are shown 
as the dashed lines in the lower panel of this figure (Kwan and Krolik 
1982). For A/B » 0.1 (eV)^ 2 our model results are similar to those of the 
standard model of KMT; for A/B = 10(eV)*/ 2 , corresponding to a B 1, our 
spectrum has roughly the same optical to X-ray flux ratio as the KMT standard 
spectrum; for this choice of spectrum we find H* about a factor of 3 larger 
than KMT's. We attribute this difference to the fact that our spectrum is 
steeper in the hard X-ray spectral region than theirs. 

The observed ensemble of radio loud quasars has a mean optical to X-ray 
flux ratio corresponding to a » 1.2 or A/B 20 (eV) *^ 2 . Radio quiet 

vA 

quasars have a » 1.5 or A/B « 110 (eV) 1 ^ 2 (Zamorani et al. 1982). Figure 1 

Ua 

shows that in either case the range H* < s < misses the inferred observed 
cloud ionization parameters entirely, although the width of the allowed 
ionization parameter region, log (-*/£*), is comparable to the observationally 
inferred width. The predicted ionization parameter range can be brought into 
accord with the observed one if the effects of cloud evaporation are 
considered. Evaporation tends to increase the cloud pressure relative to the 
ICM by a factor of ~ a few (KMT), thus shifting the clouds to lower ionization 
parameter. A shift of a factor of 6 in cloud pressure and hence ionization 
parameter, is sufficient to bring the two phase models into aqreement with the 
cloud observations. Figure 1 also shows that for A/B > 300 (eV)^ 2 , 
corresponding to a > 1.65, the two phase region is evanescent. This is 

U A 

consistent with the observed paucity of objects with spectra in this range; 

the only known measurement of a >1.9 may be a consequence of X-ray 

0 X 

absorption by intervening matter (Elvis and Fabbiano, 1984). 

Since the value of =j* is dependent on the spectrum at very high (~ 0.5 
MeV) and very low (l-10ym) photon energies, and since photons in these energy 
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ranges do not contribute directly to the cloud ionization and broad line 
emission, the correspondence between the inferred observed cloud ionization 
parameters and the computed provide constraints on the true AGN spectrum at 
very high and very low energies. Thus, our calculations rule out the 
existence of a steep (index > 1.1) IR continuum extending to wavelengths 
greater than ~ 10u in broad line objects and requires broad line objects to be 
hard X-ray sources. 

Although the two phase models predict the run of ICM temperature and 
ionization as a function of 5, the distribution of gas pressures as a function 
of position depends on a detailed dynamical model for the gas flows in AGN. 
Various suggestions for the net flow in the broad line cloud region include 
spherical accretion (e.g. Krolik and London 1982) and radiatively driven 
outflow (Mathews 1982, Beltrametti 1981 and references therein). In either 
case the gas pressure may be represented roughly as a power law function of 
distance from the central continuum source, P a R^. Free fall predicts 3 = 
-3/2, constant velocity outflow predicts 3 = -2, and radiatively driven 
outflow requires 3 < -2 (Mathews 1982). The relative merits of these models 
are influenced by the fact that the range of pressures which can explain the 
properties of observed broad emission lines do not vary by more than an order 
of magnitude over the ensemble of AGN. The range of pressures occurring 
within any given object may be much less. If we adopt the point of view that 
the broad line pressure can vary at most by a factor of 10 across the range of 
parameter space > s > s*, then 3 is determined by the relation 


3 = 


-2 

1 + Alog 5 
Alog R 


( 2 ) 
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If A log 5 = log 2* - log 2* and if 5* = 2 and 2* = 0.3 (Kwan and Krolik 
c h c h 

1982), then 3 must lie in the range -1 < 3 < 10. It is interesting to note 
that the two phase models further constrain the pressure distribution via the 
mean density of the clouds plus ICM. That is, if the mean density is such 
that the gas cannot all exist at 10 4 K throughout the region of thermal 
instability as has been shown by KMT, then the region of ionization parameter 
space 2 < 2ft must be forbidden. Thus the pressure distribution in the broad 
line region must be such that the pressure cannot be high enough for the 
ionization parameter to be less than 2ft . For a power law distribution of 
pressure with radius throughout the broad line region the power law index must 
be 3 “ -2 in order to satisfy this constraint. Finally, we note that even 
though many of these arguments concerning the broad line region density 
distribution apply to the 10 4 K clouds rather than to the ICM, the viscous 
coupling between clouds and ICM is sufficiently strong that velocity 
differences between the two components of more than ~ 100 Km sec' 1 are not 
expected. Thus these dynamical arguments apply equally to clouds and ICM, 
although hot gas may exist in regions devoid of clouds. 

Given the conflicting arguments and constraints regarding the broad line 
region pressure distribution, we attempt to explore a variety of plausible 
distributions which illustrate the range of ICM properties expected. In all 
cases we assume that the region of interest is bounded by the critical 
ionization parameters required by the two phase models, 2* and 2ft. Thus our 
estimates do not rely on the presence of hot gas in regions other than where 
it is required by the existence of clouds. 

The column density of a given ion through a spherical shell bounded by 2 
and 2 ft is 
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F J 1 


dx g(x) 

TOO 


for a power law pressure distribution, where P is the pressure at H = E£ s l is 
the source luminosity, and T is the gas temperature and g is the fractional 
ion abundance relative to hydrogen. The ion column densities are weighted 
toward ionization parameters H = if 3 > -2. Observed quasar line ratios 
indicate that clouds exist primarly near H > S* (Kwan and Krolik, 1QR2: KMT) , 
suggesting that 3 < -2 does not occur. This suggests that strongly 
accelerated winds are less likely than constant velocity (or decelerating) 
winds or virialized flow of the ICM. Since the ionization and temperature 
across the two phase region vary by a factor of only a few, the ion column 
densities for 3 < -2 are similar to those for 3 > -2. If 3 =-2, as in a 
constant velocity wind, is the true ICM pressure distribution then some other 
physical effect, such as a cutoff at large distance from the continuum source 
due to a scarcity of gas, must determine the ion column densities. Perhaps 
the simplest pressure distribution is 3 = 0, so that P = constant. In this 
case, since >> H*, the column density is: 


N « 1.5xl0 24 cm' 2 ' 


' L 46 P 14 1 

— 3* T 

*Tl '8 


(4) 


where L45 is the luminosity in units of 10 46 erg sec"*, P^ is ICM pressure in 
units of 10 14 cnf^ K, and T3 is the gas temperature in units of lO^K. 

Equation (4) provides an estimate of N for 3 not too different from zero. For 
the remainder of this section we will concentrate first on 3=0 models, and 
then argue that their results will be substantially unchanged for 3 in the 
range / 3/<2. Finally, we discuss the effects of pressure distributions which 
do not fall within this range. 


9 



c) Spectral Results 

For our choice of incident spectrum the ICM is typically highly ionized, 
and the abundance of trace ions such as the hydrogen-like species of the trace 
element 0, Si, S, and Fe depend sensitively on both the shape of the spectrum 
and on the range of ionization parameters in the ICM. The Fe XXVI and 0 VIII 
fractional abundances have maximum values, near H = Ej*, of ~ 10"* for all 
cases except for the softest spectra, A/B - 1000 (eV)^/^, for which the Fe 
XXVI and 0 VIII abundance can reach 0.4. The total column density of Fe XXVI 
is displayed in Figure 2, as a function of the spectral softness A/B through a 
constant pressure region bounded by 3* at large radius and H* at small 
radius. Also shown are the column densities of 0 VIII and the total nucleon 
column density. From this figure it is apparent that the constraint that the 
ICM be Compton thin implies that (l-4g Pi^^fA/B) 1 / 4 < 1.5. 

The cross section for absorption of photons at the center of the La 
analog resonance line of a hydrogenic ion is (e.g. Bethe and Salpeter 1957) 

(5) 

line me Ae 
e 

where e, h, m e , and c are the familiar constants, the oscillator strength in 

the non relativisitic approximation is f=0.4162, and Ae is the line width. If 

-17 2 

Ae is measured in eV then the cross section is °-|.j ne = 4.6x10 cm /A eeV. 

Thus an Fe XXVI La absorption line broadened by gas motions of 10,000 Km sec - 

*, as suggested by the UV line widths, requires an ion column density of 5 x 

10 + 1® cm"2 in order to be optically thick. The ions 0 VIII, Si XIV, and S XVI 

require column densities of 4.8 x 10^ cm^, 1.5 x 10*® cm^, and 1.9 x 10* fi 

cm^, respectively, in order to be optically thick with this choice of 
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broadening. The optical depths required for thermal ization of resonance line 
photons are much greater than those obtainable in the ICM. Therefore, if the 
gas motions are turbulent or random, absorption by material along the line of 
sight to the source will be filled in by scattered emission (the net effect of 
the loss of scattered photons by absorption in the central source is likely to 
be negligible). If the gas has a net velocity toward or away from the central 
source an observable P-Cygni, or inverse P-Cygni, line profile is expected 
(this possibility has also been suggested by Shields and McKee, 1971). The 
absorption equivalent width of such a line will be W, = /“if t a e for t, < 1 

a a a 

(Basko et al. 1980). If x a = 0.1 at Fe XXVI La we expect W g = 40 eV with our 
previous choice of broadening. Features with this strength will be observable 
with future X-ray spectrometers. Figure 2 shows that for Fe XXVI La x a =0.1 
is exceeded for all plausible choices of illuminatinq spectrum if L45 P14 = 1* 
x = 1 is attained at A/B = 30 (eV)*^. 

a 

The cross section for continuum photoabsorption by a hydrogenic ion in 
the non-relativistic limit is (e.g. Bethe and Salpeter 1957) 


0 th = 2 it 8 
cont 4 



( 6 ) 


at the threshold energy. In this equation a = is the fine structure 
constant, a Q is the radius of the first Bohr orbit, and z is the nuclear 
charge. Substituting the values of these constants gives 
°cont = 6 ’ 3xl0_18 cm 2 /z 2 . Thus column densities exceeding l.OxlO 19 cm ? , 3.1 x 
10* 9 cm^, 4.1 x lO^ 9 cm^, 1.0 x 10^ 8 cm^ are required for optically thick 
continuum absorption by the ions 0 VIII, Si XIV, S XVI, and Fe XXVI, 
respectively. Photons absorped in this way are partially thermal i zed and 
reemitted as recombination line radiation. In a recombining gas collisionally 
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excited line radiation can be neglected. The recombination line emission 


equivalent width satisfies the relation W a » t w for continuum optical 

g con c max 

depths T cont << 1. The maximum equivalent width, attained when -r cont » 1, is 


defined by 


max 


£ line ^ 


(7) 


where . is the line energy, £ is the fraction of the incident photons 
which are in the continuum of the emitting ion, and <j> is a number which 
reflects the efficiency of conversion from photoabsorption to line emission. 
For our range of spectra 0.1 < £ < 0.4, and <j> ~ 0.2 for case A recombination 
emission of La {e.g. Pengelly 1972). Thus, t c > 0.05 is required in order 
to produce W e > 10"^ e ij ne * From figure 2 it is apparent that Fe XXVI is 
(again) the best candidate for an observable feature. T con t = 0.1 is exceeded 
for A/B > 50 (eV)*/2. Thus, any recombination emission line will be weaker 
than an absorption feature or any scattered emission in a P Cygni profile by a 
factor ~ 10. However, recombination line emission will appear even if the qas 
flow is predominantly turbulent or random. 

We now explore the consequences of possible pressure variations across 
the two-phase regions. Figure 3 shows the result of evaluating equation (3) 
using ion fractions and temperatures as a function of E calculated for an A/B 
= 10(eV)l/2 spectrum. This figure shows that the conclusions reached so far 
are not qualitatively changed when pressure variations are considered. For 3 
> - 2 the column density is influenced predominantly by conditions at H = 
and increments of 3 result in only small diminutions of the column 
densities. For example, Np eXXyi never decreases below 10 17,7 (L^ P^r^cm -2 
for - 2 < 3 < + 2, corresponding to a La analog line optical depth of at least 
0.1 using our previous broadening estimate. 
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HL CONSEQUENCES OF TWO PHASE MODELS: COMPTON THICK CASE 

We now consider the case when the ICM is optically thick to Compton 
scattering. In the context of the two phase models for the broad line region, 
as outlined in § II, this can arise if the gas pressure, source luminosity, or 
source shape is such that the column density through a uniform sphere bounded 
by H = 2* (c.f . equation (1)) exceeds 1.5 x 10 24 cm -2 . Our treatment of this 
problem is influenced by the fact that the gas pressure may be greater closer 
to the source than in the broad line region. In addition, the two phase 
models provide only a lower limit for the ICM density--the ICM may be more 
dense than required by pressure equilibrium (see e.g. Eilek and Caroff 
1982). Thus, rather than restricting our choice of ICM properties based on 
the two phase models, in this section we examine the general effects of 
Compton scattering by gas near a power law continuum source. In § IV we 
discuss a possible scenario whereby this can occur. The effects of Compton 
scattering on the broad line profiles have been discussed by Shields and McKee 
(1981); we will ignore them here. 

The primary effect of multiple Compton scattering on an AGN power law is 
a reduction in the flux at the highest photon energies, relative to the rest 
of the spectrum. The mean fractional energy shift per scattering of a photon 
of energy e is (e.g. Ross, Weaver, and McCray, 1978). 


Ae> _ 4kT e 
e ~ me 2 ” me 2 


( 8 ) 


in the non relativistic limit. For a gas temperature T ~ 10^ K approximately 
150 scatterings are needed to increase the photon energy appreciably, while 
only ~ 2 scatterings are needed to significantly decrease the energy of 250 
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keV photons. We will restrict ourselves to optical depths and gas 
temperatures such that the effects of photon energy gain due to Compton 
scattering can be neglected. 

The consequences of Compton softening for the thermal equilibrium in the 

broad line region is illustrated in Figure 4. This figure shows the locus of 

equilibrium temperatures for a gas illuminated by a spectrum which has passed 

through a uniform (density) sphere of Compton depths 0, 4, and 8. The initial 

spectrum used here is our standard double power law (Equation 1) with A/B = 10 

(ev)l/2; the Compton softening was calculated using the procedure described by 

Illarionov, et al. (1979); the thermal equilibrium was calculated in the same 

way as in the previous section. For t c = 4, the softening causes the spectrum 

to have an X-ray to optical ratio resembling that of the softest observed 

AGN. For t > 4 the spectrum is softened sufficiently to entirely suppress 
C + 

v 

the thermal instability. Objects with Compton depths in excess of this value 
may not have broad line clouds, unless an additional heating mechanism 
exists. Gas near these objects will be distributed uniformly at a temperature 
of 10 5 - 10 6 K. 

A second consequence of multiple Compton scattering is that the photon 
path length in the scattering medium will be much greater than the size of the 
medium. This increases the probability of photoabsorption if the medium is 
not fully ionized. In the 10 5 K - 10 6 K gas at the boundary of a Compton 
thick cloud the photoelectric opacity will be much greater than if the qas 
were Compton thin and the temperature were 10 7 - 10 8 K. The amplification of 
this absorption by multiple Compton scatterings can suppress the escaping flux 
even if the opacity is such that the gas would be optically thin to absorption 
in the absence of multiple scatterings. Furthermore, the steeply decreasing 
AGN power law spectrum will preferentially photoionize the low energy 
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absorbers H, He, and C. Thus the gas can be transparent between 13.6 eV and 
the K edge of 0 VIII, near 800 eV (Krolik and Kallman 1984). 


A summary description of the effects of simultaneous Comptonization and 

photoabsorption by ICM gas is as follows; (1) Photons lose a fraction of their 
Ag 

energy <— > at each scattering. The total energy loss on traversing a medium 

o Ag 

of depth t c is approximately <— > (e.g. Illarionov, et al . 1979). (2) The 

total photon path length in the medium is approximately t c R, where R is the 

size of the medium. The optical depth to photoabsorption is therefore x g = t c 

R K a , where * a is the absorption opacity. (3) The local photon flux in the 

medium is increased by a factor proportional to x, owing to the increased path 

length. In order to model these effects with high accuracy it is necessary to 

simultaneously solve the equation of radiation transfer, including the effects 

of energy redistribution due to Compton scattering, together with the 

ionization and thermal equlibrium equations. Such a procedure is cumbersome 

and not well suited to probing a wide parameter space {e.g. Ross, Weaver, and 

McCray 1978). We therefore adopt an approximate technique for modeling 

simultaneous Comptonization and photoabsorption which preserves all of the 

qualitative behavior of more accurate treatments without necessitating a 

detailed solution of the full radiative transfer problem. We do so by 

modifying the program used to calculate the thermal equilibrium results 

described in §11 in the following ways: (1) At each radius step in our 

spherical clouds a fraction of photons at energy e are removed, the fraction 
Ac 2 2 

being <— :> dx where dx is the increment of Compton depth squared in that 
radius step. (2) The local ionizing flux in the interior of the cloud is 
increased over the Compton thin value by a factor 1 + x c , where x c is the 
Compton depth to the cloud boundary. (3) The photoabsorption opacities are 
increased over the Compton thin values by a factor of 1 + x c . As a further 
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computational expedient, the clouds are assumed to have constant density 
rather than pressure. This allows the effect of the cloud boundary to be 
calculated analytically in the cloud interior, and avoids the need for an 
iterative solution of the transfer and thermal equilibrium. The appropriate 
ionization parameter in this case is proportional to the ratio of the 
radiation flux F at the cloud boundary to the gas density n, which we define 
as 5 = 4irF/n (Tarter, Tucker, and Salpeter, 1969). The properties of such a 
Compton thick cloud are determined entirely by the Compton depth x c> 5, and 
the shape of the incident continuum. 

We have tested our approximate treatment of Comptonization by comparing 
the spectrum escaping a pure scattering cloud of depth t c = 4 calculated with 
our technique with that calculated using the formulation of Illarionov et al . 
(1979). The approximate spectrum is found to be within 20% of the exact 
result for photons which have been shifted by more than ~ 0.5 Compton 
wavelengths as they escape. For small wavelength shifts, corresponding to 
photon energies greater than ~ 100 KeV, our approximate technique 
significantly underestimates the escaping flux. We attribute this to 
predominantly non-Gaussian nature of the distribution of energy shifts for one 
or two scatterings. The error involved is small since even the exact 
formulation predicts that less than 10% of the flux emitted above ~ 100 keV 
escapes in that energy range from a Comptonizing cloud of t c 2. 

A sample Comptonized and photoabsorbed spectrum is shown in Figure 5, for 
x c = 10 and log 5 - 2.5. The incident continuum is shown as the dashed curve, 
and the effects of simultaneous Comptonization and photoabsorption are shown 
in the solid curves. The effects of photon destruction at the highest 
energies by Compton scattering and in the 1-10 keV band by photoabsorption are 
apparent in this figure. Photons below ~ 800 eV escape the cloud freely for 
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this choice of parameters. The energy removed from the X-ray band above - 1 
KeV escapes the scattering and absorbing region at lower energies, either as a 
result of down-Compton scattering or as reprocessed photo-absorbed radiation 
(i.e. as soft X-ray and XUV lines and continua from ions such as 0 VIII, C VI 
and He II). The detailed dependence of the escaping spectrum on cloud depth 
and the ionization parameters at the cloud boundary is displayed in Figure 
6. In this figure the contours of constant flux in the 13.6 eV - 500 eV 
energy band (dashed lines) and 500 eV- 4 keV band (solid line) are shown as a 
function of t c and 5. These fluxes are measured relative to their optically 
thin values, for a spectrum of the form given by equation (1) with A/B = 10 
(eV) 1 / 2 . 

For 1 < ? c < 2 the shape of the escaping spectrum is a sensitive function 
of t Cj owing to the degredation of the highest energy photons (e - me") by 
Comptonization. Since the gas temperature is coupled to the high energy flux 
via Compton heating, the ionization state and opacity are also sensitive to t c 
in this region. If we define 6 C as the critical ionization parameter below 
which photoabsorption becomes important, the value of K c (inferred from the 
constant flux surfaces labeled 1.0 in Figure 5) is sensitive to x c , for 
x c < 2. For 2 < t c < 10 the effects of Comptonization are complete: further 
increases of x c have negligible effect on the escaping spectrum and on 6 C . 

The effects of preferential ionization of low energy absorbers in this range 
of x c is apparent from the fact that the 0.5 - 4keV flux decreases to 30% of 
its optically thin value for by log (5) < 2.2, while the 13.6 eV - 0.5 keV 
flux is less strongly affected unless by log (5) < 1.8. The effects of 
Compton amplification of photoabsorption are apparent for t c > 10. In this 
range of depths the photons in the 0.5-4 keV range can be 90% absorbed, 
while less than 50% of the 13.6 eV - 0.5 keV photons are absorbed. This can 
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occur over a range of ionization parameter of roughly a factor of 10 at t c = 
20 as illustrated by the hatching in Figure 6 . This behavior is consistent 
with the results of analytic estimates for the transmission coefficient of a 
Compton scattering, photoabsorbing sphere. Such estimates (see appendix) 
predict that surfaces of constant transmitted flux will be congruent 
with t a £ surfaces, and that photoabsorption will be more important at high 

v» 

photon energies than at lower photon energies. 


IV. DISCUSSION 

In addition to accounting for the existence of quasar broad line clouds, 
the two phase models can be used to predict ICM properties for AGN as a 
function of continuum luminosity. The cloud filling factor, f, is expressed 
in terms of cloud density, n c , and total (cloud + ICM) mass M accordinq to: 

(9) 

T" R um H "c 

where v and m^ are the gas mean molecular weight and proton mass, 
respectively. The broad line region size is given by the ionization parameter 


R = 


f L il/2 
''4ttcP = ' J * 


( 10 ) 


so the filling factor is 


f - 1.7 X 10’ 3 L-J M 3/3 P 14 =. ( 11 ) 

where L 45 is the source luminosity, M 4 is the mass in units of 10^ Mg and P^ 
is the pressure in units of 10 14 cm - 3 K. Thus the filling factor will approach 
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unity for L ~ 10^3 erg sec - *, in accord with the observed almost total 
covering of low luminosity Seyfert galaxies (Holt, et al . , 1980; Reichert et 
al. 1984), and the low mean covering factor in OSO's (Osmer 1980). 

The combined column density of the clouds and ICM is given by M = n R, 
where n, the mean gas density, is given by (KMT) 

n = n c f + (1-f) n h (12) 

and n c and n^ are the cloud and ICM gas densities, respectively. The column 
density is therefore 


N = 


3M Pc = 


n. 


u m, 


H 




1/2 


kT 


IC 


(13) 


if the pressure is constant across the broad line region. This equation shows 
that the column density increases proportional to (LP) 1//2 at high luminosity, 
in accord with equation (4). However, the two phase model also predicts a 
column density increase proportional to (L/P ) “^ at low luminosity, if M, and = 
are independent of luminosity. Thus, we expect Compton thick ICM gas both at 
very high luminosity, and at very low luminosity. This is consistent with the 
idea that Seyfert 2 nuclei are lower luminosity versions of Seyfert 1 1 s , and 
that their continuum and broad line properties are dominated by the effects of 
Comptonization and Compton-amplified photoabsorption. This is also consistent 
with the apparent cutoff in the quasar luminosity function above - 10^ erg 
sec -1 (e.g. Kuhr et al . 1983). Very high luminosity quasars are likely to 
have their emission lines smeared into unobservability by Compton scattering. 
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as noted by Shields and McKee (1981), in addition to having a weakened or 
absent thermal instability. 

The observed properties of Seyfert 2 galaxies include weakened or absent 
broad emission lines (Weedman 1976), weak or absent X-ray flux in the 0.5-5 
keV ( Einstein Observatory) X-ray band (Lawrence and Elvis 1982), and weak or 
absent X-ray flux at energies greater than 5 keV (Piccinotti et al. 1982). 
These galaxies also show narrow emission lines which resemble those of higher 
luminosity Seyfert 1 galaxies. The Comptonized and photoabsorbed spectra 
presented in § III have sufficient flux escaping in the 13.6 eV - 800 eV 
energy range to penetrate a column density ~ 10^ cm"^ in narrow line clouds 
and produce narrow lines similar to those observed (Halpern and Steiner, 

1983). 

Other observed properties of Seyfert 2's support the idea that they are 
optically thick to Compton scattering. Optical polarization studies 
(Antonucci, 1983) have revealed that the linear polarization of Seyfert 2's is 
preferentially perpendicular to the radio structure in these objects, while 
Seyfert 1 1 s show polarization which is preferentially parallel. If the gas in 
these galaxies is distributed in a flattened cloud whose normal is coincident 
with the radio structures, then the observed polarizations are consistent with 
optically thick and optically thin Compton scattering for Seyfert 2's and 
Seyfert 1 ' s , respectively. Although our discussion in this paper is based on 
an assumed spherical gas distribution, the qualitative results will be 
unchanged if the distribution is flattened slightly. As shown by Angel 
(1969), an axial ratio of 1/3 is adequate to account for the degree of 
polarization of a few percent observed from Seyfert 2's. Finally, we note 
that the optical continuum of Seyfert 2's exhibits variability limits that are 
less than those of Seyfert 1 ' s . This behavior is consistent with the effects 
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of Compton scattering smearing variability on timescales less than the broad 
line region light crossing time, ~ 1 week. 

Our model results show that a Compton depth of ~ 1 is sufficient to 
suppress the two phase instability, while a depth of 5-10 is required to 
suppress the escaping X-ray flux at > 1 keV. Therefore, if the ensemble of 
known active galaxies has a distribution of Compton depths spanning this 
range, the range of AGN properties should include those with both broad lines 
and simple power law X-ray spectra (e.g. quasars and Seyfert l's), those with 
neither broad lines or power law X-rays (e.g. Seyfert 2 1 s ) , and those with 
weak or absent broad lines but apparently normal X-ray spectra. Possible 
examples of this last class of objects include the intermediate Seyfert's, and 
the narrow emission line galaxies (Mushotzky, 1982). The observed 
metamorphosis of some intermediate type Seyfert's into objects with no broad 
lines may also be consistent with the effects of Comptonizaton. For example, 
the dynamical time scale at the radius of the broad line region in the Seyfert 
galaxy NGC 4151 is roughly 1 year. On this timescale an influx of material to 
the broad line region could increase the Compton depth sufficiently to 
suppress the thermal instability and account for the observed disappearance of 
the broad lines (Ulrich, et al. 1984). It should be noted that the X-ray 
spectrum and the UV/X-ray flux ratio are not observed to vary stronqly in NGC 
4151 (Perola et al. 1984). 

The results of this paper can be summarized as follows: In §11 we showed 

that two phase models predict ICM column densities in excess of 10^ cm - ^ for 
a wide range of plausible AGN conditions. This conclusion depends only on the 
requirement that the observed broad line clouds be confined by ICM thermal qas 
pressure. Furthermore, the degree of ionization of the ICM is predicted to be 
sufficient to produce potentially observable X-ray line absorption. For 
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example an Fe XXVI La absorption equivalent width of 40eV will occur for 
almost any plausible illuminating spectrum and ICM pressure distribution if 
the ICM pressure and source luminosity satisfy the condition > 1. 

Features of this strength will be observable with future X-ray spectroscopy 
experiments, such as those planned for the Advance X-ray Astrophysics Facility 
(AXAF) satellite. The absorption lines will only occur if the ICM gas has a 
net velocity toward or away from the continuum source. Emission lines will 
have equivalent widths approximately a factor of 10 smaller than for 
absorption lines, and will be produced even if the gas flow is turbulent or 
random. Detection of X-ray lines from AGN provides a test for the two phase 
models ofthe broad line clouds, and also serves as a probe of the broad line 
gas motions. 

In §1 1 1 we showed that Compton scattering by the ICM can have two effects 
which will affect both the broad line and X-ray properties of AGN. At 
moderate Compton optical depths the spectral softening due to Comptonization 
can suppress the two phase instability which allows the coexistence of broad 
line clouds and ICM in quasars and type 1 Seyferts. At larger Compton optical 
depths (t > 10) scattering increases the photon path length sufficiently to 
cause absorption of photons in the 1-10 keV energy band. Such a medium will 
also be transparent to photons with energies less than 1 keV if the ICM mass 
is comparable to that inferred for quasars. We therefore suggest that this 
scenario accounts for the absence of X-rays and broad lines in Seyfert 2 
galaxies. This idea is consistent with the polarization and variability 
properties of Seyfert 2's as compared to Seyfert l's. It also arises 
naturally if the cloud plus ICM mass and ionization parameter are assumed to 
be the same for all AGN. 
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APPENDIX 


Analytic Theory of Photoabsorption with Compton Amplification 

The equation of ionization equilibrium for an ion whose ionization is 
controlled primarily by its own photoabsorption is: 


, _L_ f de f -x{ £ ) 

( o - n 2 J e. u e T e e 
4ttR th 


°( £ ) = X + n a (T) 


(Al) 


where x 0 is the fractional abundance of the ion in question, x+ is the 
franctional abundance of the next higher ionization stage, L is the total 
incident luminosity, R is the distance from the source, f e is the spectral 

oo 

distribution function (/ f £ d £ = 1), a is the photoionization cross section, 

n e is the electron number density, and a(T) is the recombination rate 

coefficient. If we assume that: (i) x+ " 1, Xq << 1: (ii) T = const; 

(ill) ! ^ f £ e _T(e) a(e) - f e e' T{e) O(e) £=£th then 

th 


4ira x 
7 al e 


(A2) 


or 


e _T dx 


4nan 

t r 


dR 


(A3) 


if n = constant this integrates to 


i q - t _ 4 Tra (Ln) 1//2 
1 " e ~ eT“ ~^J2 


(A4) 


The transmission coefficient of a Compton thin cloud of Compton depth x c and 
ionization parameter C at the boundary is therefore 


24 




(A5) 


T = 


-T 

e 


T 

4ita C 


where is the Thomson cross section. In the Compton thick case the 
effective photoabsorption cross section is inreased by a factor t c , so the 
transmission coefficient becomes 


T = 


e 


-TT 

C 


4ira T C 


(A6) 


The important conclusions which can be reached from this discussion are 
that: (i) Surfaces of constant T are congruent to surfaces of t c a C, i.e. 

diagonal lines, in the t 5 plane; (ii) ions with smaller specific flux in 
their ionizing continuum, i.e. ions with higher threshold energies, will be 
more likely to recombine and cause photoabsorption than will ions subject to 
larger ionizing fluxes. Thus at a given value of t c and 5, ions with high 
energy photionization threshold (e.g. e > 1 keV) are more likely to cause 
optically thick photoelectric absorption than will ions with low energy 
threshold (e.g. £ < 1 keV). 
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Figure 1(a) Compton Temperature (T IC ) as a function of the spectral softness parameter, A/B as 
defined m equation (1) 

Figure 1(b) - The limiting ionization parameters bounding the two phase region, and hJJ, plotted 
as functions of the spectral softness parameter, A/B Dashed lines denote the 
allowed region of ionization parameter space as inferred from emission line ratios 
and photoionization models for the broad line clouds (KMT). The value of A/B 
which corresponds to the maximum value of a ox is shown as a dashed line 
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Figure 6. Contours of constant flux in the 13 6 eV - 0 5 keV (dashed) and 0 5 keV-4keV 
(solid) energy bands as a function of cloud Thomson depth, t c , and boundary 
ionization parameter, £. Curves are labeled with the fluxes relative to the incident 
spectrum, which is assumed to be the same as in Figure 5 In the hatched region the 
emergent spectrum will resemble that shown in Figure 5 
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